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Abstract 

When poured into a container, cohesive granular materials form 
low-density, open granular aggregates. If pressed upon with a 
ram, these aggregates densify by particle rearrangement. Here 
we introduce experimental evidence to the effect that particle re- 
arrangement is a spatially heterogeneous phenomenon, which oc- 
curs in the form of a phase transformation between two configu- 
rational phases of the granular aggregate. We then show that the 
energy landscape associated with particle rearrangement is consis- 
tent with our interpretation of the experimental results. Besides 
affording insight into the physics of the granular state, our con- 
clusions are relevant to many engineering processes and natural 
phenomena. 

Cohesive granular materials have been the focus of a minute fraction of recent 
research into the granular state. 1 Yet cohesive granular materials will surely draw 
increasing attention from scientists and engineers, if only because they are used 
in numerous applications. Conspicuous examples are afforded by the forming of 
ceramic parts, pulvimetallurgic components and pharmaceutical tablets by com- 
paction of fine powders. 2 The cohesiveness of powders stems from the large sur- 
face/volume ratio of their constitutive particles, which enhances the effect of at- 
tractive van der Walls forces. In other applications, e.g. the stabilization of soils, 
the cohesiveness is due to the presence of liquid menisci among the particles. Our 
interest in the densification of cohesive granular materials was prompted by the 
recent compaction study of Kong and Lannutti. 3 These authors used X-ray to- 
mography to document the evolution of density during the compaction of alumina 
powders (diameter ~ 60 /mi). They reached the tantalizing conclusion that den- 
sification "seems to proceed as a wave initiated at the advancing ram". 3 Our aim 



Figure 1: a) Computer simulation of container filling by a ballistic aggregation method. 8 
The particles are sequentially deposited along random vertical paths, and then allowed to 
roll until they make contact with any two points (at least) on the surface of the aggregate. 
The magnitude of the cohesive forces is assumed strong enough to stabilize a particle 
both at the time of aggregation and afterwards, but it remains otherwise unspecified, b) 
Obtained particle aggregate, c) More or less regular rings of particles are a pervasive 
feature of the simulated aggregate, and also of d) aggregates obtained in quasi two- 
dimensional container-filling experiments. Note the water menisci among the particles, 
e) Collapse of a ring of particles by snap-through buckling: one of the particles jumps 
to the center of the ring. The driving force F is a relatively large contact force localized 
on the jumping particle; both photo-elastic and numerical studies indicate 5 that such 
forces exist with magnitudes several times larger than the average contact force, which 
is represented by the hydrostatic pressure p 



here is to elucidate the nature of this 'wave,' and to relate its behavior to the 
micromechanics of densification in cohesive granular materials. 

When, preceding compaction, a cohesive granular material is poured into a 
container, the mobility of the particles reaching the bottom of the container is 
hindered by the cohesive forces, Fig. |l|a. As a result, a low-density, open aggre- 
gate of particles obtains inside the container, Fig. [I|b. Open aggregates densify 
by particle rearrangement at relatively low pressure. 4 It has been proposed 5 that 
particle rearrangement occurs when the rings of particles of the open aggregate 
collapse by snap-through buckling, see Fig. |l]c-e. To investigate this phenomenon 
we prepared a quasi two-dimensional open aggregate 6 by filling a narrow plexiglass 
container (thickness ~ 1.9mm) with monosized glass beads (diameter ~ 1.7mm). 
Before pouring them into the container, we wetted the beads with water in order 
for menisci to form among the beads. These menisci provided the required cohe- 
sion. We then compacted the aggregate using a ram. Fig. ^|a shows three stages 
during the experiment. A high-density region (the phase H, wherein rearrange- 
ment has taken place already) and a low-density region (the phase L, wherein the 
open aggregate 

remains unchanged) are clearly discernable. Densification proceeds by growth 
of the volume fraction of H at the expense of the volume fraction of L. Visual 
inspection revealed that no rearrangement occurs within the high- or low-density 
regions. In a narrow vicinity of the H-L interface or rearrangement front, on the 
other hand, we could clearly observe the collapse by snap-through buckling of suc- 
cessive layers of rings of particles, see Fig. ^|b. The rearrangement front, which we 
identify with the densification wave of Kong and Lannutti, advances in the form 
of broad, shallow tongues (visible in Fig. ^|a) darting forward in coincidence with 
the collapse of ring of particles. These results suggest that densification occurs in 
the form of a phase transformation L^H, Fig. |^a-b. To substantiate out interpre- 
tation of the experimental evidence we turn our attention to the micromechanics 
of particle rearrangement. Consider a ring of particles undergoing snap-through 
buckling, Fig. |I]e. For an increasing displacement A, Fig. ^c, the force F vanishes 
as the particle snaps into the void, and again when the particle reaches the center 
of the void. Fig. |3|d shows the attendant evolution of internal energy. The energy 
associated with the initial (tangent) response of the ring, W t , is represented by a 
dashed curve; the relaxation effected by buckling, W^, is represented by an arrow. 
Wb causes the function W(A) to be nonconvex. Nonconvex energy functions are 
characteristic of systems which undergo phase transformations. 7 

To move on to the macroscopic scale we consider an open granular aggregate 
contained in a frictionless container of volume V and constant cross-sectional area. 
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Figure 2: a) Three sucessive times during densification. The ram is black, and moves 
from the top down. The high- and low-density regions are dark and light gray, respec- 
tively. The area shown is ~ 20 x 20 cm 2 . b) Clockwise starting from the top left: 
successive collapse of rings of particles in a vicinity of the rearrangement front. The 
particle diameter is 1.7 mm. 
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We study the energetics of densification in the space of the local stretch A. (The 
local stretch is defined by A = po/p, where p is the local density and po is the 
density of the initial open aggregate.) The energy per unit volume is 0(A), Fig. [3]e; 
it is nonconvex by inheritance from W(A), Fig. ||d. The local pressure is given by 
p = —d(j)(\)/d\. We now set the average stretch to a given value, A = y J v XdV, 
and minimize the total energy of the aggregate, J v <f)(\)dV, using conventional 
tools of nonconvex analysis. 7 This leads to the following equilibrium equations 



which allow for the computation of the characteristic stretches Al and Ah and the 
transformation pressure p t , Fig. ^f. It is apparent from Eq. [I] that Al, Ah and p t are 
independent of A, and can be construed as material properties. The characteristic 
stretches define two configurational phases L and H of density p L = Po/Al and 
Ph = Po/Ah, respectively. When Al > A > Ah, the phase H occupies a volume 
aV, and the phase L a volume (1 — a)V, where a = (Al — A)/ (Al — Ah), Fig. |]a. 
As a increases from to 1 during densification, A decreases fom Al to Ah, the 
average density increases from pl to pn, and the pressure remains equal to p t . It 
is noteworthy that p t ^ on account of the dissipative effects of inter-particle 
friction, Fig. |3|e-f. The average energy per unit volume is given by the convexified 
form = «0(Ah) + (1 — «)0(Al), which for Al > A > Ah yields the correct value 
p = —d(f){X) I d\ = p t . 

The spatial distribution of phases is not given by our analysis. Nuclei of H could 
form in many places throughout the aggregate, leading to stratified mixtures of H 
and L. There is no indication of such mixtures in Fig. |]a, however. We ascribe 
this fact to the regularizing effect of container- wall roughness. Because of wall 
roughness, the stress decreases monotonically away from the ram. Nucleation 
must therefore occur at the ram, leading to the establishment of a single H-L 
interface — the rearrangement front. We will explore other effects of wall roughness 
in a separate publication. 

We conclude that if, as indicated by our experimental observations, snap- 
through buckling is the dominant mechanism of particle rearrangement, then den- 
sification must occur in the form of a phase transformation L^H. Container-wall 
roughness accounts for the presence of a single H-L interface. It follows that our 
interpretation of the experimental evidence is consistent with the micromechanics 
of particle rearrangement. 
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Figure 3: a) Interpretation of the experimental results: V is the total volume of the 
aggregate; the phases H and L are separated by the rearrangement front; a is the volume 
fraction of H. As a increases from to 1, the rearrangement front sweeps through the 
aggregate from ram to bottom. Cf. Fig. ||a. b) The density jumps by Ap = pn — Ph 
across the rearragement front, c) Mechanical response of the ring of particles of Fig. |l]e: 
force F vs. displacement A (schematic), d is the particle diameter, d) Internal energy 
W vs. A. e) We express the energy per unit volume of aggregate as <j) = (f>t + (f>b + 4>f , 
where <f>t and <f>i correspond to the terms Wt and W& of (d) averaged over a statistically 
representative volume of aggregate, and 0/ is the energy per unit volume dissipated 
by the inter-particle frictional forces, f) Geometrical interpretation of the equilibrium 
equations, Eq. |l[ g) The average energy per unit volume 4>(X). 
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